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ABSTRACT 

Forty-five high-latitude, OB stars have been observed In the Lya and 21 cm 
lines of HI in an effort to map out the vertical distribution and extent of 
the local HI halo. The 25 stars for which a reliable HI column density ^an be 
obtained from Lya lie between 60 and 3100 pc from the plane. Our principal 
result is that the total column density of HI at |z| >1 kpc is, on the 
average, 5 ± 3 x 10*9 cm“2, or 15% of the total At relatively low |z| 

the data toward some stars suggest a low effective scale height and fairly 
high average foreground density, while toward others the effective scale 
height is large and the average density is low. This can be understood as the 
result of irregularities in the interstellar medium. A model with half of the 
HI mass in clouds having radii of a few pc and a Gaussian vertical 
distribution with a z » 135 pc, and half of the mass in an exponential 
component with a scale height of 500 pc, gives a satisfactory fit to the data. 
The technique of comparing Lya and 21 cm column densities is also used to 
discuss the problem of estimating the distance to several possibly subluminous 
stars. 

Subject headings : galaxies: Milky Way - galaxies: structure - 

radio sources: 21 cm radiation - stars: luminosities - 

ultraviolet: spectra 
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I. INTRODUCTION 

The combination of Lya absorption measurements, which give the amount of 
interstellar neutral hydrogen toward a star, and 21 cm emission measurements, 
which give the HI toward and beyond a star, can be used to map out the 
distribution of neutral gas in the galactic halo [fiohlin, Savage and Drake 
1978; Hobbs et al . 1982 (Paper I)]. The main advantage of this technique is 
its independence of a priori assumptions about the form of the vertical 
distribution of hydrogen. The main difficulties are practical: the lack of 

suitable high latitude target stars for Lya measurements, observational 
uncertainties in N a , and the Inaccuracies in 21 cm column density 
determinations caused by stray radiation, beam dilution and opacity 
corrections. In Paper I we examined data toward nine distant, high-latitude 
stars of early spectral type, and concluded that the neutral halo extended at 
least 1 kpc from the plane and possibly further. In this paper we reconsider 
the question of the vertical extent of halo HI using a larger group of target 
stars and a new set of sensitive 21 cm observations that are free of stray 


radiation. 
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XI. OBSERVATIONS 

(a) The Program Stars 

The A5 stars for which we have Lya and 21 cm data are listed in Table 1. 
The stars were chosen to satisfy the following requirements: (1) |b| > 30° , 

to minimize the effects of gas in the galactic disk; (2) an MK. temperature 
class of B2.5 or earlier, to minimize the effects of the stellar photospheric 
contribution to the Lya line; (3) apparently accurate spectral types, to 
provide reliable stellar distances and to implement, reliably, the temperature 
discrimination intended in item 2; and (A) 6 > -A5° , to allow suitable 21 cm 
observations from Green Bank. Table 1 includes nearly all stars at 6 > -25° 
which meet these criteria and have V < 6.5 and therefore contains a 
considerable fraction of all such stars with M/M 0 > 6, owing to the small 
scale heights of 0 and B stars. 

The first two requirements were later relaxed slightly. Our sample 
includes two stars at b = -27°, p Col, and HD 20366A, from the classic survey 
of Munch and Zirin (1961). Twenty of the stars, furthermore, have spectral 
types between B2.5 and B5 Inclusive, and indeed eventually proved of little 
use for our purposes. They were originally included on the chance that their 
nominal spectral types were incorrectly too late and, In some cases, because 
earlier studies of their interstellar line spectra had been reported. 

The color excesses E(B-V) in Table 1 are derived from the intrinsic colors 
of Johnson (1963) and the observed spectra and colors given in the table. ^Jew 
photoelectric B, V magnitudes and spectral types obtained at Yerkes 
Observatory were kindly made available. In advance of publication, by 
K. M. Cudworth and W. W. Morgan, respectively. Values of E(B-V) for Be stars 
are marked with a colon, since the emission lines may affect the observed 
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colors and thereby introduce relatively large errors into the derived color 
excesses. The spectroscopic distances, d, in Table 1 are obtained from Lesh's 
(1968) calibration of MK spectral types and include corrections for 
interstellar absorption calculated by assuming Ay/Eg-y «* 3.1. The distances 
differ somewhat from those of Paper I owing to the choice here of Lesh's 
(1968) absolute magnitude scale in place of Blaauw's (1963). The values of d 
and z marked with a colon are more uncertain than the others owing to larger 
uncertainties in the spectral types: either these spectra are peculiar, or 

the types were derived from high-dispersion spectra not intended and not 
optimally suitable for MK classification. Radial velocities with respect to 
the LSR, V r , listed in the next to last column of Table 1, are derived in most 
cases from the heliocentric values of Hoffleit (1982) or of Munch and Zirin 
(1961). 

(b) The Lya Observations 

The column densities N a (HI) obtained from various Lya observations are 
given in column 6 of Table 2. Twelve of the values are previously published 
results from Copernicus observations. The remaining values are taken from the 
recent IUE survey of Shull and Van Steenberg (1985; hereafter SVS). Most of 
these stars were incorporated into that survey specifically for the present 
investigation, and the results were obtained in about equal numbers from new 
IUE observations and archival IUE spectra. The basic method of analysis in 
all cases was that of line-profile fitting to the radiative-damping wings of 
the observed Lya line, as described, for example, by Bohlin, Savage and Drake 
(1978) or by SVS. A typical observational uncertainty in the final value of 
N a is ± 0.1 dex, according to both sets of authors. We note that SVS have 
observed 11 of the 12 stars in Table 2 for which we have listed the N a values 
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from Copernicus . For most of these stars the N a (HI) value obtained from XUE 
agrees with the Copernicus value to within the stated errors; in all cases the 
agreement is within twice the stated errors. 

A colon follows entries in columns 6 and 8 of Table 2 if the star is of 
spectral type later than B2. The strong stellar Lyct line present in these 
cooler stars makes a reliable estimate of the Interstellar HI column density 
generally impossible (Savage and Panek 1974; SVS). In most cases, this effect 
makes the apparent values of N a /E(B-V) for the later type stars considerably 
greater than the "standard’' N(H)/E(B-V) ■ 58 x 10^0 cm“2 mag"^ found by 
Bohlin, Savage and Drake (1978) and SVS. Stars of spectral type B2 or 
earlier, however, have values of N a /E(B-V) in agreement with 
58 x 1020 cm-2 mag -1 within the uncertainties in N a and a plausible 
uncertainty in E(B-V) of 0.03. For these reasons, analyses carried out below 
use Lya data from only the 25 stars of spectral type B2 or earlier. Any 
spurious increase in the value of N a by the stellar Lya line in all but a few 
of these stars is expected to be much smaller than a factor of two. 

(c) The 21 cm Observations 

The 21 cm data were obtained during 1983-84 at the telescopes of the NRAO 
in Green Bank. The primary instrument was the 140' (43 m) telescope which has 
a half-power beam-width (HPBW) of 21' at 21 cm. It was used to measure the HI 
profile toward each star at least three independent times in the course of 
four observing sessions., There were two different receiving systems. One had 
scalar feeds and a system temperature T S y S =* 50 K. The other had hybrid mode 
feeds and T S y 6 = 23 K. Most stars were observed with both systems. The 
velocity resolution was always 1 km s"*, the spectra covered ± 250 km s~^ 
about zero velocity LSR, and the average rms noise in the final spectra is 
0.06 K (T b ). 
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All of the radio spectra initially contain stray radiation, that is, 21 cm 
emission that comes from arbitrary parts of the sky but is scattered or 
diffracted into the receiver. When the telescope is pointed at high galactic 
latitudes, where the hydrogen column density is small, stray radiation can 
constitute more than half of the observed signal (see, e.g., van Woerden, 
Takakubo and Braes 1*362; Giovanelli et al . 1978; Kalberla, Mebold and Reich 
1980). Stray radiation was removed from the spectra using a method described 
by Lockman, Jahoda and McCammon (1985). In brief, an area ~10 square degrees 
around each star is mapped in HI with the 140' telescope. The maps are then 
convolved to the '"3° x 2 ° angular resolution of the Bell Labs HI survey (Stark 
et al . 1985) . Because the Bell Labs survey is virtually free of stray 
radiation, differences between the convolved 140’ data and the Bell data are 
interpreted as stray radiation which is then removed from the Individual 140* 
spectra. Numerous tests indicate that this technique gives reliable, 
repeatable estimates of the HI brightness temperature. 

The amount of stray radiation that was removed from 140’ spectra varied 
with the stellar coordinates and the hour angle of observation, but it was 
sometimes equivalent to column density corrections of -1 x 10 20 cm -2 and thus 
makes a significant difference in the estimated column density N 21 toward some 
stars. For example, we find only 76% and 55% of the column density reported 
by Giovanelli et al . (1978) toward p Col and HD 93521, respectively, and 89% 
of the column density used in Paper I for HD 93521. 

The main source of uncertainty in the 21 cm spectra, after removal of 
stray radiation, is instrumental baseline error. This was estimated 
individually for each star from the differences between the 21 cm profile for 
the various observing sessions. The total range of fluctuations Is given as 
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the error estimate in Table 2. It is always larger than the errors due to 
noise in the data. 

There is a final uncertainty in the conversion of 21 cm brightness 
temperature to HI column density because this requires knowledge or the spin 
temperature of each component in the HI spectrum. It is not possible to know 
the spin temperature without making emission/absorption measurements (see, 
e.g., Dickey, Salpeter and Terzian 1978; Payne, Salpeter and Terzian 1981; and 
Dickey and Benson 1982), and these require that the direction of interest 
coincide with a strong extragalactic radio source. Table 2 gives column 
densities for two assumptions about the hydrogen: that it is optically thin 

(infinite spin temperature), or that it has a T s = 75 K. In many directions 
the choice makes little difference to N2J,. We will use the T a ® 75 K values 
for comparison with N a because they are, in the median, 5 % larger than the 
optically thin values, and this is comparable to the Increases obtained from 
studies at high latitude (Dickey and Benson 1982). The emission/absorption 
experiments cited above indicate that opacity corrections are always ,£10% at 
| b | > 25°. 

Some of the stars were also observed with the 300' (91 ra) transit 
telescope at Green Bank at an angular resolution of 11’ and a velocity 
resolution of 1.4 km s“l. Because it is awkward to map a large declination 
range with the 300 r telescope, these data cannot be corrected for stray 
radiation by direct comparison with the Bell Labs HI survey. Instead, when 
the 300' data are used (only to investigate possible beam dilution effects due 
to the relatively large solid angle covered by the 140' observations), a stray 
radiation correction was estimated by forcing the 300' column densities, 
averaged over the area of the main 140* beam, to agree with the corrected 140' 
column densities. 
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III. RESULTS 

(a) Are the Directions of the Stars Typical of the ISM? 

Before making a comparison of Lya and 21 cm data it is important to 
determine if there is anything out of the ordinary about the Interstellar gas 
in the directions of the distant, high latitude, early-type stars. For this 
purpose the program stars are divided into two groups. The "primary sample" 
consists of the 25 stars with spectral type earlier than B2.5 that are 
presumed to have reliable Lya column densities. The "control" stars are those 
with a later spectral type and generally unreliable estimates of N a . The 
control group also includes 21 cm data toward 12 stars, mostly of spectral 
type A, that we observed in the course of other work, for a total of 32 
control stars. A comparison, unbiased sample of interstellar HI is given by 
the Arecibo 21 cm data in the direction of ~50 extragalactic radio sources 
(Dickey, Salpeter and Terzian 1978; Payne, Salpeter and Terzian 1982; Liszt. 
1983). The main quantity of interest is the vertical component of the HI 
column density, which, for the three samples, is 

<N 2 i sin|b|> =* 3.2 - 1.9 x 10 20 cm -2 (extragalactic) 

<N 2 i sln|b|> a 3.3 ± 1.4 x 10 20 cm -2 (primary sample) 

<N 2 i sin | b | > => 3.4 I 2.3 x lO 2 ^ cm” 2 (control stars) 

where all the uncertainties are la. These numbers show that there is nothing 
at all peculiar about the total Ngj in the direction of either set of stars; 
they seem typical of the ISM. We note, parenthetically, that there are no 
high-velocity clouds (|v| > 100 km s“*) In the 21 cm spectra toward the 
primary sample stars, although a few such clouds appear within several degrees 


of some of the stars. 



The agreement between the average of the Green Bank data, which were 
obtained at 21' resolution, and the average of the Arecibo data, which were 
obtained at an angular resolution ~4’ indicates that the larger beam gives a 
faithful sample of <N 2 i> and its dispersion. In addition,” the 21 cm HI maps 
that cover several square degrees around each star (SXIc) allow us to compare 
N 21 in the direction of a star with N 21 in nearby regions of the sky, and to 
examine beam-to-beam fluctuations in N 2 i* Each star has a measured N 21 C*) 
and, in a circle of 66' radius centered on the star, a mean <N 2 i> and its 
dispersion, o. Averaged over each group, 

N 2l(*) / <N21> " 1 • 0 1 ± 0.06 (primary sample) 

N2i(*) / <N 21 > “ 1.00 ± 0.07 (control stars) 

and 

a/< N 2 i> M 0.12 ± 0.04 (primary sample) 
o/<N 21 > " 0.12 ± 0.05 (control stars). 

There is neither a surplus nor a deficit of HI toward the stars compared to 
their immediate neighborhood, and the column density fluctuations in that 
neighborhood are comparable to the HI column density fluctuations found on 
angular scales from tens of arc-sec to a few degrees (Jahoda et al . 1985). 
Also, there is no significant difference in the 21 cm emission between the 
direction of the target and control stars, despite the difference in spectral 
type of the groups. In sum, there is no evidence in the 21 cm data that the 
younger stars are located in peculiar surroundings or that their directions 
are atypical of the interstellar medium at high latitudes. Furthermore, there 
appear to be negligible amounts of molecular hydrogen at high latitudes, at 
least along the path to the 10 of our program stars studied by Savage et al . 
(1977), so the HI results should pertain to the entire neutral medium. 
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b) Comparison of Lya and 21 cm column densities 

The data give three quantities of interest.' the observed ratio N 2 i/N a , 
the observed difference (N 2 i*-N a )sin|b | , and an average "foreground" density 
<n> ■ N a /d e ff, where d R ff is an effective path length of gas. Somewhat 
arbitrarily we adopt d e ff (z) - h[l-exp(-|zj/h)]/sin|bj where h is a 
characteristic thickness of the HI layer. Thus the effective distance is in 
the range 0 < d e £f < h/sin |b| for all z. With this definition, the average 
density <n> is identical to the central density n(0) in the Idealized case of 
a uniform exponential layer With a scale-height of h (Paper I). Average 
densities are usually calculated from <n> * N a /d, i.e., with d e ff * d, but 
this choice is less satisfactory for our sample of distant, high latitude 
stars, some of which probably lie well beyond most of the HI layer (at 
z » h) . Use of d e ff * d for these stars can give values of <n> that depend 
more on the star's distance than on the average HI density in the disk. 

The observed quantities are shown in Figures 1 and 2, where the ratio and 
the difference are plotted against |z|. Error bars are ± la and are dominated 
by uncertainties in N a . The ratio of 21 cm to Lya column densities vs. |z| is 
drawn again in Figure 3, with two additional kinds of information. First, ' 
each star is represented by its value of <n>, in units of 0.01 cm - ^, 
calculated for the illustrative value h ■ 500 pc. Second, the curve shows the 
ratio N 2 i/N a for a uniform, exponential disk with a scale height of 500 pc. 
Points that lie above the curve indicate equivalent exponential scale heights 
h e q >500 pc, while those below the curve imply values h R q <500 pc, where 
h e q (z,r) = |z| {£n[r/(r-l) ] and r = N 2 i,/N a . In contrast to the quantity 
<n>, which is independent of N 21 and therefore characterizes only the gas in 
front of the star, the quantity h e q depends on the ratio N 2 j./N ct at any z. It 
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therefore concisely combines N^i/Nq and z to characterize the diotribution of 
the gas both in front of and behind the star. 

The main conclusions from these figures arc: (1) The uncertainties in the 

»• 

data are so large that the amount of high-z HI is not well determined. The 
column density beyond |zj > 1000 pc could be negligible, or perhaps as much as 
' 1 x 1020 cm“2, (2) Stars at similar |z| do not always have identical 

values of the ratio or difference. (3) There are significant differences in 
<n> between stars, even if they are at the same |z[ (e.g., the stars around 
|z| - 300 pc in Fig. 3). The difference can be more than an order of 
magnitude — far larger than the experimental uncertainty. There is more 
dispersion in <n> at low |z| than at high jz|. (4) There is a trend, 
especially at low |z|, for stars with high average density <n> to have low 
ratios of N 2 i/N a , indicating a small equivalent scale-height for the gas, 
while low densities are associated with a large equivalent scale-height, 

(5) There is one star, HD 214930, that appears to have N a > N 21 at a 
significant level. This is most simply explained as stellar contamination of 
the interstellar Lya line, and indeed, the star has a N a /E(B-V) that is twice 
the usual value. But a N a > N 21 can also occur if a very small, very dense 
cloud lies toward the star and suffers severe beam dilution in the 21 cm 
measurements. 

The remainder of this section is devoted to the five main points derived 
from Figures 1-3. These suggest that the data cannot be Interpreted in terms 
of a smooth HI layer, but demand a multi-component medium with different scale 
heights for each component. 
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c) Angular Resolution Effects and Structure In the ISM 

It Is posolble that the vast difference between the volume of the 
Interstellar medium (ISM) sampled by a 21 cm emission measurement and that 
sampled by a Lya absorption measurement is the source of much of the structure 

in Figures 1-3. There is an extensive body of 21 cm data and analysis that 

addresses the question of small scale fluctuations. It can bo summarized in 
the statement that the total column density of high-latitudc HI 1 b remarkably 
smooth on small angular scales. Specifically, there is no significant 
difference between the statistical properties of HI seen against pulsars, 
arc-sec sized extragalactic radio sources, or extended extragalactic radio 
sources. There are not large differences in the total Ngi toward the separate 
components of extragalactic double radio sources. At high latitudes the total 
HI column density is highly correlated on scales ranging from < 30" to several 
degrees, and the values of o/ <%x> ~ 0.1 derived for the vicinity of the 
program stars is typical of radio observations in general. Of course, the ISM 

is not uniform, and some structure is observed on all angular scales. But 

fluctuations in total column density on scales < 2° are modest, usually only 
"■ 10% of the mean. The idea that a significant fraction of HI is contained in 
sub-parsec size features (the "cloudlet" or "subcloudlet" model), producing 
an extremely mottled HI surface density across the sky, is now thoroughly 
discredited. Significant changes in across the sky result from gradients, 
not from small-scale random fluctuations (Dickey, Salpeter and Terzian 1978; 
Dickey 1979; Dickey et al . 1981; Dickey and Benson 1982; Liszt, Greisen and 
Dickey 1982; Payne, Salpeter and Terzian 1983; Jahoda et al . 1985; Dickey, 
Crovisier and Kazes 1984) . In view of this there is every reason to expect 
that a star well above the HI layer will have N a = N 21 with a ler uncertainty 


10 % 
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The data in Figure l show that N a is within 10% of N 21 (± la) at 
|z| >1 kpc. The weighted average of (N 2 i-N a )sin|b | for the six stars at 
| z | >1 kpc is 0.5 i 0.3 x 10^0 cm" 2 . Thus, within the experimental error, 
the data are consistent with an HI layer that does not extend substantially 
beyond 1 kpc, and is well behaved in that fluctuations are equal to those 
expected from a priori considerations. A reasonable estimate is that 
NHlsinjb) < 1 x 1020 cm" 2 at |z| > 1 kpc. This is the principal result of the 
present work. 

In contrast to the most distant stars, stars at low jz| should show wide 
variations in the ratio N 2 i/N a and in <n> because about one-half of the HI 
mass is contained in "clouds" (e.g., Falgarone and Lequeux 1973; Spitzer 1978; 
Liszt 1983). Two stars embedded in the HI layer at Identical |z| could have 
very different values of N a , and therefore of <n> and h e q, if one lay 
behind a cloud (or, for that matter, behind a hole) and the other did not. 

This appears to be the case for stars at |z| < 800 pc, where the fluctuations 
show up as discordant values of N 2 j/Na and <n> for different stars at the same 
1*1 . 

In summary, the Lyot and 21 cm data are consistent with previous knowledge 
about the general features of the HI layer: it has a vertical extent not much 

greater than 1 kpc, and there are many irregularities in its lower regions, 
d) Model Calculations 

Calculations that simulate the Lya-21 cm comparison can clarify the 
conclusions of the previous section on the effects of a multi-component medium 
and the difference between the angular resolution of the Lya and 21 cm 
measurements. The model discussed here gives the HI column density toward 
test stars located at the distance and latitude of the observed stars, and 
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also gives the beam-averaged N 21 In those directions. The radio observations 
were simulated by determining column densities at '*400 positions around each 
star and averaging them with a weight appropriate to the radio telescope beam. 

The model contains the minimum ISM necessary to reproduce the data: a 
smooth, low density HI component and a set of denser clouds. Because the 
equivalent <n> toward the stars is ,£1 cm“^ we did not Include a separate class 
of "large" clouds in the simulation (see Spitzer 1985) . The clouds we use are 
homogeneous and spherical with r'.tner a constant radius or a radius drawn from 
a power law distribution; N(r c ) * r c "2. All clouds In a given simulation 
have the same HI volume density. To keep the model as close as possible to 
known properties of interstellar HI, the clouds contain 50% of the mass, they 
are located randomly but have a Gaussian vertical distribution with a 
dispersion, o z , of 135 pc, and the number of clouds was adjusted so that the 
average line of sight would intersect 8 kpc”^ in the plane (e.g., Falgarone 
and Lequeux 1973; Crovisier 1978; Liszt 1983, and references therein). No 
clouds were allowed within 20 pc of the sun, although the results are not very 
sensitive to the radius of this local hole. 

The "smooth" component is less constrained by existing observations; if 
does not even have to be smooth. The main requirement is that it have a very 
much higher spatial frequency than the clouds. Since there are, on average, 
eight clouds kpc" 1 in the plane, and the half-thickness of the HI layer Is 
only a few hundred pcs, many high latitude directions would have no HI if 
there were nothing but clouds. This is not observed. The smooth component 
was approximated by an exponential layer with a scale height of 500 pc (after 
Bohlin, Savage and Drake 1978; Lockman 1984). When given a central density of 
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0.1 era" 3 (e.g., Spitzer 1985), this component contains half of the observed 
<N 2 i ein|b|> ~ 3 x 10 20 cm“ 2 , and is consistent with our data (§IIIc) for 
| z| >1 kpc. 

The only free parameter in the model is the cloud size. Although many 
simulations were run with a power-law distribution of radii, the uncertainties 
in the data make only the grossest features of the models testable. For this 
reason we will discuss only models in which all clouds have the same radius: 
0.5 < r c < 20 pc. 

Results from several hundred simulations show that clouds with radii of a 
few pc make the HI column density fairly well correlated on angular scales 
< 2° but completely decorrelated on larger scales, where the fluctuations 
converge to o/<N 2 ]> ~ 0.3. Averaged over the sky this type of model will 
always be somewhat smoother than the 21 cm observations mainly because it has 
no large-scale gradients. On angular scales ~20* , the model has fluctuations 
o/<N 2 i > ~0.36 r c _1 . Clouds with radii > 2.5 pc give a plausible 
approximation to the irregularities in N 21 and also are the appropriate size 
for the smallest major fluctuations (Dickey, Crovisier and Kazes 1984), so 
the discussion will be restricted to these. 

Figure 4 shows N 2 i/N a when the model ISM, with r c = 2.5 pc (and thus 
n c « 12.5 cm -3 ) is observed with an infinitesimally small antenna beam 
(HPBW = 0) . The points include both representative and extreme values from 40 
simulations. The curve shows N 2 i/N a for the 500 pc scale height of the 
exponential component — the thickest component in the model. Points above the 
curve imply equivalent scale heights h e q > 500 pc, those below the curve imply 
h e q < 500 pc. 
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Many features of the simulation correspond to the observations. First, 
there Is a large range in the ratio N 2 i/N a for test stars at low z: the Btars 

at 200 < | z | < 300 pc, for example, have ratios between 1.2 and 4.2. The 

t * 

largest ratio comes from a simulation which produces a single cloud .lust 
beyond the star; the smallest ratio comes from a simulation in which 6 clouds 
lie toward the star but none beyond it. The variety occurs despite the 
identical angular resolution of the radio and UV simulations. Second, because 
the column density through a cloud is comparable to the column density through 
the entire exponential layer, many stars at |z| < 1 kpc show the dominance of 
the clouds to the exclusion of the exponential component. This shows up as a 
high <n> coupled with a low N 2 i/N a . These model stars (typified in our 
observations, we think, by 53 Ari and £ Lib) give low values for the effective 
layer thickness while overestimating its average density. Third, at 
|z| > 300 pc the simulated ratio never lies above the curve marking 
h e q = 500 pc. The model cloud layer has a vertical dispersion of 135 pc, so 
stars at |z| > 300 pc are beyond almost all of the clouds. Fourth, at large 
|z] the ratio approaches unity and the density is in the range 0. 1-0.4 cm - ^. 
Over large distances both the 21 cm and Lycc data faithfully reflect the 
average values of the model and the dispersion about the average caused by the 
clouds. Finally, because there is no beam dilution in this simulation, N 21 is 
always greater than N a , and thus the ratio is never < 1. 

Figure 5 shows the same simulations as Figure 4 when N 21 is obtained at a 
HPBW of 21’ . The individual values of N a remain unchanged. The only obvious 
effect of the large beam is an increase in the scatter of the ratio, 
especially at large z. There are several ratios < 1 that occur when clouds 
obscure the star but do not fill the radio beam, and the ratio can also lie 
above the h e q = 500 pc curve at large |z|. Apart from the increased scatter, 



I 


18 

differences between the HPBW ■ 0 and HPBW ■ 21* models are, in general, much 
smaller than the uncertainties in the data. The mean absolute difference 
between N21 at HPBW - O' and HFBW - 21' for this model is only 

1 » 

0.2 ± 0.2 x 10 20 cm -2 , while the maximum absolute difference is 
0.8 x 1020 cm" 2 . 

Most features of the observations are reproduced in Figure 5. The model 
gives a range of the ratio and <n> at low z, gives a general anticorrelation 
between the ratio and <n>, and converges to a ratio of unity and an <n> of a 
few tenths at large |z|. This model also is in reasonable agreement with the 
observed differences (Fig. 2). 

A number of simulations were run with different values of the exponential 
scale height, h, and the cloud layer thickness, c z , to see if the observations 
could constrain these parameters. The results are rather imprecise owing to 
the limited amount of data and the relatively large errors. However, the 
exponential component fits the data best if it has 500 < h < 750 pc. With 
lower values the models hardly ever give ratios like those observed toward 
y Peg, p Col and especially p Leo 5 with greater values there are too many 
stars at low z with very large ratios. The model with h = 500 pc gives a 
%isin|b| of 2 x 10 19 cm -2 for |z| > 1 kpc, in agreement with our 
measurements. For the clouds, a layer with o z = 200 pc does not place enough 
gas at low |z| to give the high densities observed toward 53 Ari and 5 Lib, 
while with a z - 70 pc the layer is so thin that there are only a few dense 
clouds at high Latitude; stars tend either to have a ratio -1 or to lie on the 
h e q = 500 line. 

In conclusion, a very simple two component model is in reasonable 
agreement with most aspects of the data. The "standard” model with a Gaussian 
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cloud distribution of a z - 135 pc and an exponential "Intercloud" distribution 
with h - 500 pc, gives a satisfactory fit. This model is derived from, and 
hence in good agreement with, other observations, 
e) KD 214930 

In the direction of HD 214930, which has a spectral classification of B2 
IV (Table 1), the ratio N 2 i/N a - 0.55 ± 0.10 (la). This low ratio Is most 
simply explained as stellar contamination of the Lya line, as would be 
expected if the true spectral type were somewhat later than B2. This 
explanation would also account for the relatively high value of N a /E(B-V) 
toward the star (Table 2). But, if the low ratio is not the result of stellar 
Lya absorption, or of an undetected experimental error, then it can only be 
caused by the presence of a dense cloud that obscures the star but is of such 
small angular size that it does not contribute much to N 21 . The properties of 
this cloud are firmly constrained by the observations. 

Toward HD 214930 we have 21 cm data at 21* and 11’ resolution. The two 
21 cm HI profiles are identical to within the noise: the higher resolution 

spectrum has the same N 21 as the lower resolution spectrum and shows no new 
line components. If there Is a small HI cloud toward the star then its column 
density, averaged over the 11' beam, must be < 2 x 10^ cm -3 (the magnitude of 
the uncertainty in N 21 ), while the column density through the cloud in the 
direction of the star must be ~4.5 x 10^0 cnf 3 to account for the difference 
between N 21 and N a . If the hypothetical cloud is homogeneous, spherical, 
and directly in line with the star, and if there is no HI at all beyond the 
star, then the cloud, located |z c | pc from the plane, must have a radius r c 
and a density n c of 

r c (pc) < | z c | / 1200 ; 
n c (cm~ 3 ) > 9 x l£)V|z c |. 
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The value of |z c | is limited to < 520 pc by the star; therefore r c < 0.4 pc, 
and n c > 170 cnf^. It must have a diameter < 3*. These limits are firm, 
because any change in the assumptions used to derive them will shrink the 

i • 

radius further and raise the density unless the cloud is extended exactly 
along the line of sight. 

The HI in the field around this star is fairly smooth. Over an area of 
50' x 40' fully sampled at 5' intervals with an 11' HPBW, the column density 
fluctuations o/<N 2 i> are only 0.08, and much of the dispersion is caused by a 
gradient of *>1 x IQ 20 cm" 2 per degree that runs from west to east. The column 
density of 21 cm HI toward the star is virtually identical to the average over 
the field: N 2 i(*)/<N 2 i> ■ 1.03. Thus there is no indication of oignificant 
angular structure in N 21 over this area of the sky. 

The cloud that would be required to reconcile the 21 cm and Lya data seems 
quite improbable. First, it is so dense that it should be predominantly 
molecular rather than atomic. Molecular gas is rare at high latitudes (Savage 
st al » 1977; Blitz, Magnani and Mundy 1984) and, when present, is typically in 
clouds '■*1° across and not in very small clouds. Second, the "small cloud" 
hypothesis implies fluctuations in N^x of 100% over angular scales of a few 
arc-min. As stressed in §IIIc, fluctuations of this magnitude have not been 
observed. These considerations, combined with the fact that HD 214930 appears 
to have a gas-to-dust ratio about twice the usual value, suggests that Na for 
this star is too high, most likely because of stellar contamination of the 
interstellar line component. 
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IV. POSSIBLE SUBLUMINOUS STARS 

In Paper I we gave the detailed reasons for our conclusion that the 9 
stars under consideration were true, massive, Population I objects. In view 
of those reasons, our current primary sample of relatively bright stars (which 
Includes the original 9) should contain few, if any, subluminous stars of the 
type that are spectroscopically indistinguishable, at classification 
dispersions, from unevolved, massive Pop I stars of similar temperature. It 
then follows that the spectroscopically determined distances d and z in 
Table 1, and our present conclusions based on them, are generally reliable. 
However, the incidence of low-luminosity stars in any sample of much fainter 
high-latitude blue stars is an important and incompletely understood issue. 

A valuable ll3t of six candidate subluminous stars with <V> 9.6 has 

recently been reported by Dworetsky, Whitelock, and Carnochan (1982; hereafter 
DWC) in their Table 7, which specifically excludes the larger numbers of 
spectroscopically recognizable sdO and sdB stars also found by those authors. 
These six candidate "old disk subdwarfs" appear to be anomalously unreduened 
at the distances inferred from their apparently normal spectra, according to 
DWC. Any independent means of estimating the distances to these stars, and 
hence their absolute magnitudes, is important in assessing this conclusion. 
These six stars were therefore considered for inclusion in our observing list, 
on the possibility that their distances might be constrained from their ratios 
N2i/N a . Four of the six candidates, however, were not suitable for 
observation because of their relatively late spectral types, low latitudes, or 
excessively long estimated IUE exposure times. They were replaced with six 
other candidate disk subdwarfs, taken from Table 8 of DWC, that have no 
available classification spectra, but typically appear, from their colors, to 
be slightly hotter than the original candidates. 
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The observational resultB are collected in Table 3, where the stellar data 

in the first six columns are from DWC. The values of N a in column 8 have also 

been reported by SVS. In contrast to our primary sample of stars (Table 2), 

* «* 

for which 96% have a ratio N 2 i/N a > 0.9, only 2 of 8 (25%) of the DWC stars do 
so, despite equal observational uncertainties for the two groups. This 
statistically well-marked difference indicates that the stellar contributions 
to the Lya lines of the DWC candidates in Table 3 typically exceed the 
interstellar contributions, and therefore that most of the stars are somewhat 
cooler than inferred by DWC. This conclusion is strengthened by the fact that 
the two stars for which MK types were determined directly by DWC, HD 121800 
and HD 148265, do have ratios N 2 l/N a consistent with their spectral types. 

No firm, new, conclusions as to the possibly subluminous nature of the 
stars in Table 3 can be drawn from the N 2 i/N a ratios alone. As an extreme 
example, if a number of the stars actually are B5V Population I objects with 
My « -1 and B-V = -0.16, then typically d ~ 1500 pc, |z| ~ 800 pc, 

E(B-V) ^ 0.03, and their main-sequence lifetimes of ~6 x 10? yr are 
sufficiently long to allow them to have been born in the plane with low 
initial velocity (|v z | < 20 km s”l). As emphasized by DWC and in Paper I, an 
independent astrometric estimate of these distances, in the form of very 
accurate proper motions, is needed. 
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V. DISCUSSION 

The data in Paper I suggested that as much as one-third of the HI at high 

latitudes was between 1 and 2 kpc from the galactic plane. The new data 

»# 

presented here, improved both in quantity and quality over what was previously 
available, modifies this picture of the halo considerably: the total column 

density above |z| “1 kpc now seems to be only 15 % of the total, and may be an 
even smaller fraction. 

The new, expanded data set also shows the complex vertical structure of 
HI, Below | z | ~ 800 pc there are directions where the equivalent HI scale 
height appears to be quite low, and also directions where the scale height is 
> 500 pc. There are also variations in the average density, <n>, by a factor 
"•40 among stars at low jz|, but only a factor ~4 among stars at high |z[ 

(Fig. 3). Density fluctuations do appear to average out over large distances, 
but at low [ z | a large number of stars must be observed to obtain a complete 
picture of the medium. 

The data suggest an ISM that is, in spirit though not necessarily in 
detail, like that in the two component model of cloud and intercloud (Clark 
1965; Field, Goldsmith and Habing 1969; a good recent discussion of this is in 
Liszt 1983). Indeed, a model with a Gaussian cloud layer and exponential 
intercloud medium gives a reasonable, though hardly unique, fit to the data. 
Because the program stars were selected to be at high latitude, we naturally 
avoid lines of sight through the denser parts of the ISM, and our data do not 
sample the full range of interstellar densities. 

Up to this point we have not discussed our observations in the context of 
the McKee and Ostriker (1977) model of the interstellar medium. This model 
envisions a neutral ISM that is entirely in clouds with dense cores and 



tenuous outer envelopes. The overlapping envelopes mimic an intercloud 
medium. Our data arc not well fit by such a model. First, there is so much 
hydrogen in the small cores that there are large fluctuations in the total Ngx 
on small angular scales. Such a medium would give a much larger scatter in 
N2l/ N a> even l z l >1 kpc, than is observed (see also §IXIc). Second, the 
observations of £ Lib and 53 Arl, which show a low effective scale height and 
high average density for the gas, require that substantial amounts of HI be 
located at low |z|. But a population of core-halo clouds that meets this need 
will not simultaneously give the large effective scale height and low <n> 
observed toward p Col and p Leo without violating (at the minimum) the 
observed limits on N 2 isin|b|. Our data are most easily fit with two 
populations that have different scale heights and average density. This 
criticism of the McKee-Ostriker model has been voiced previously by Helles 
(1980) and Liszt (1983) in their discussion of the distribution of HI as 
derived from radio observations. While it may be possible to construct an ISM 
entirely of discrete entities whose properties vary with |z| so ac to match 
all observational requirements, and there is even a suggestion that HI clouds 
change with |z| in the necessary sense (Crovisler 1978), a model of this sort 
has not yet been presented. 

To summarize, (1) Most of the HI is below l kpc, though there can be 
significant amounts above 500 pc, as seen, for example, in the line of sight 
toward p Leo. This is consistent with the general picture of the vertical 
distribution of HI that has emerged from UV studies of the local region 
(Bohlin, Savage, and Drake 1978; SVS) and radio studies of the inner Galaxy 
(Lockman 1984). (2) There is good agreement, over long paths, between Ngi and 
N a , indicating that the fundamental technique of Lya-21 cm comparison is sound 
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despite the difference in angular resolution between the observations, and 
also that small-angular scale fluctuations in the total %j are modest. (3) 
Irregularities in the HI layer at jz| < 500 pc are so large that many stars 

t « 

must be observed to define the characteristics of the layer. Our current 
sample is not adequate for this purpose. More target stars in the range 
400 < | z | < 1000 pc would be especially valuable, for at these distances there 
should be sufficient HI beyond the stars to be easily detectable, while the 
density fluctuations in the lower layers of the ISM should have averaged out. 


We acknowledge financial support fotf part of this research from the 
National Aeronautics and Space Administration, through grant NAG 5-255 to two 
of us (L.M.H. and J.M.S.), as well as contributions by C. E. Albert and 
M. Van Steenberg to portions of the data reduction and analysis. The IUE data 
analysis was supported by NASA grant NAS5-26409. 
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APPENDIX A 
Tha 21 cm spectra 

Many of our program stars have a long history as targets for ISM studies 
(e.g., Munch and Zlrtn 1961; Albert 1983), The 21 cm spectra that wo 
obtained, being of high sensitivity and essentially free of stray radiation, 
may be useful in future work. They are given in Figure 6. Spectra in the 
21 cm line toward any of the stars listed in Table 1 may be obtained from the 


authors. 
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TABLE 3 

Candidate Disk Subdwarf b 


HD/BD 

SL 

b 

V 

B-V 

MK 

*21 
» • 

T s = 75 K 

(1 0 20 cjij-2) 

Ncx 

(1020 cm “2) 

N 2 i/N a 

+30°0057 

117° 

-31° 

9.60 

-0.11 

• • • 

6.7 

8.9 

0.75 

+32°0270 

133 

-29 

10.29 

-0.16 

• 1 • 

4.5 

6.5 

0.69 

+19°0302 

142 

-40 

9.99 

-0.10 

• II 

7.5 

10. 

0.75 

121800 

113 

+50 

9.11 

-0.17 

B1.5V 

1.9 

1.5 

i.3 

148265 

45 

+42 

9.68 

-0.11 

B2.5VS 

4.1 

7.6 

0.54 

+ 7°4795 

69 

-37 

10.67 

-0.12 

• • • 

6.1 

11. 

0.55 

+16°4689 
„ a 

77 

-31 

10.37 

-0.15 

• • • 

4.6 

8.9 

0.52 

+31 4958 

107 

-28 

10.87 

-0.13 

• • a 

5.5 

5.5 

1.0 


a 2 
The blue star labeled HOB 1 by SVS, which lies about /'NW of BD +3T4958. 
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FIGURE CAPTIONS 

Figure 1. The ratio N 2 i/N a vs |z| for stars with apparently reliable Lya 
column densities. Error bars are ± la. 


Figure 2. The difference between the 21 cm and Lya column densities corrected 
by sin|b|. Error bars are ± la. 


Figure 3. The curve shows the ratio N 2 i/N a expected from a uniform 
exponential layer with a scale height of 500 pc. The numbers give the average 
density, <n>, a3 defined in the text, in 0.01 crn*^, for the observations. 


Figure 4. The ratio N 2 i/N a vs. |z| for a simulated two-component interstellar 
medium observed with an infinitesimally small antenna beam (HPBW » 0). The 
smooth curve shows the ratio if the HI were in a uniform exponential layer 
with a scale height of 500 pc. Points are shaded according to their average 
density, <n>, derived in the same way as for the observations. The test stars 
displayed are the extreme and representative values from 40 simulations. 


Figure 5. The same model simulations as in Figure 4, except that the 21 cm 
data were calculated for a 21' HPBW radio beam. 


Figure 6. The 21 cm spectra toward the program stars. 
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